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intermediate gauge group between the standard model and the scale of unification of the three coupling con- 
stants. We shall assume that there is no necessarily a group of grand unification at the scale of convergence 
of the couplings. With this scenario, different 331 models and their corresponding supersymmetric versions are 
considered, and we find the versions that allow the symmetry breaking described above. Besides, the allowed 
^ ' interval for the 331 symmetry breaking scale, and the behavior of the running coupling constants are obtained. 

, It worths saying that some of the supersymmetric scenarios could be natural frameworks for split supersymme- 

try. Finally, we look for possible 331 models with a simple group at the grand unification scale, that could fit 
' the symmetry breaking scheme described above, 
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1 Introduction 

Since the birth of the Standard Model (SM) many attempts have been done to go beyond it, and solve some of 
the problems of the model such as the charge quantization and the unification of the gauge couplings. In some 
cases the unification is done by taking a simple group of grand unification, arising the so called Grand Unification 
Theories (GUT), where the three interactions described by SM are treated as only one ^ ne mos t common 

GUT's are 50(10) and E$. The first condition for these kind of theories is an equal value for the three couplings at 
certain scale of energy, Mjj. This condition cannot be fulfilled by the simplest grand unification schemes with the 
minimal SM particle content and taking the precision low-energy data. However, the minimal supersymmetric SM 
can achieve this scenario for the coupling constants @]-[B]. Other possibilities for unification are the introduction 
of more degrees of freedom like fermions and scalar fields that lead the three couplings to converge at a high energy 
scale [TJ. Polychromatic extensions of the SM i.e. SU(N)c <8 SU(2)l <g) U{l)xm have been considered where 
the unification of gauge couplings is achieved with N — 7, 5 and with two, three Higgs doublets, respectively 
[S]. Alternative proposals of unification of quarks and leptons at TeV scale were considered too [5]. Finally, 
another interesting alternative consists of enlarging the electroweak sector of the SM gauge group such that the 
renormalization group equations (RGE) could lead to unification of the three gauge couplings at certain scale 
Mjj, in which there is no necessarily a group of grand unification at the scale of convergence of the couplings. In 
particular, the model based on the SU(3)c <£> 5(7(3) l®U {l)x gauge group (hereafter 331 models) is an interesting 
choice that could address problems like the charge quantization [TD]-[TT] and the existence of three families based 
on cancellation of anomalies [T5]- 117| . 

One way to look for new Physics in 331 models is to check for tT production with t denoting the ordinary 
top quark and T being an exotic quark with charge 2/3 or 4/3 according to the model considered. In LHC the 
production channels would be pp — > X° — » tT, when the charge of T is 2/3, and pp — > X ++ — > tT, in the case 
in which T has an exotic charge 4/3. In the first case, in order to identify the signals, the decay of t is known 
and identified by the energy spectrum and angular distribution of the final fermions |18j-[20j. then this signal is 
correlated with the one of an exotic particle with the same charge of the top, but with a totally different decay in 
the final state. The decays of T would be of the form T — > X°t — > vEt, T — > K + b — > vEb, which could de easily 
identifyable if exotic leptons have already been produced at LHC. In the second case in which T has an exotic 
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charge 4/3, it should be looked for decays of the type T — > X ++ t where X ++ is a doubly charged gauge field that 
should be easily identified if this model is correct. The other possible channel is T — > K + b — > vEb. 

In the scenario of SUSY 331 models, different channels can be searched. A good perspective is pp — > gg — > g — > 
5^ where 3 — ► iT, tT, TT in the case in which T posseses an ordinary charge 2/3. Another production mechanism 
in the gaugino sector is could be pp — > A++ — * x ++ X°: ~~ * X ++ ~~ * A++x°, — * Z, Z' ++ \ [2"T]-[2"5] 

In looking for unification of the coupling constant by passing through a 331 model, we shall assume that 1) The 
331 gauge group is the only extension of the SM before the unification of the running coupling constants. 2) The 
hypercharge associated with the 331 gauge group is adequately normalized such that the three gauge couplings 
unify at certain scale Mjj. and 3) There is no necessarily a unified gauge group at the scale of convergence of the 
couplings Mjj . In the absence of a grand unified group, there are no restriction on My coming from proton decajQ. 

Under our scheme, we have three characteristic energy scales: Mjj where the three gauge couplings converge, 
Mx where the 331 symmetry is broken, and Mz where the SM breaking occurs. We are going to consider different 
scenarios for 331 models with one and three families^], i. e., one family models. We also introduce Supersymmetric 
versions of the 331 models with different scalar Higgs multiplets. As for the SUSY breaking, we shall consider two 
different scenarios: when SUSY is broken at the electroweak scale, or when SUSY breaks at the scale of Mx- 

In our scheme we have four parameters to take into account, and to look for a possible unification of the coupling 
constants (UCC). They are the scales Mx, My, the value of the coupling contants at the unification convergence 
point, ay, and the parameter associated with the normalization of the hypercharge (denoted by a). Since we 
are interested mostly in possible phenomenological scenarios, the relevant parameter will be the gauge symmetry 
breaking scale Mx; and the parameters My, a can be viewed as functions of this one. 

If the unification came from a grand unified symmetry group G, the normalization of the hypercharge Y would 
be determined by the group structure. However, under our assumptions, this normalization factor is free and the 
problem could be addressed the opposite way, since the values obtained for a could in turn suggest possible groups 
of grand unification in which the 331 group is embedded, we shall explore this possibility as well. 

In the present work, we study six different versions of 331 models with non-SUSY and SUSY particle content, 
and find which models could lead to a unification at certain scale My with only one symmetry breaking between Mz 
and Mjj scales. Some of the SUSY versions studied, could provide a quite natural scenario for split supersymmetry. 
Finally, we also consider the possibility of embedding those 331 versions into a grand unified theory (GUT) in which 
a gauge group at the unification scale appears. 



2 Running Coupling Constants 

The evolution for the Running Coupling Constants (RCC) at onedoop order is ruled by the solution of the Renor- 
malization Group Equations (RGE), which can be written in the form [28J: 

1 (2.D 



where ctfj = gf/4w, and the coefficients 6j are given by [2T)j : 

h = + lzZ T ^ s ) - y^(G). (2.2) 

/ s 

The summations run over Weyl fcrmions and scalars, respectively. The coefficient Tr is the Dynkin index 

Tr{T a T b ) R = T R 6 ab , (2.3) 

with the generators in the representation R. The last term is the quadratic Casimir for the adjoint representation 

f acd f bcd = C 2 (G)6 ab , (2.4) 

with 02(G) = N for SU(N). On the other hand, the respective supersymmetric versions are 

b susY = Tm (f) + J2 Tm(s) - 3C 2l (G). (2.5) 
/ 

where the usual non-supersymmetric degrees of freedom are counted. 



1 Notwithstanding, proton decay could be induced even in the absence of a group of grand unification. It may occurs via six 
dimensional 331 invariant effective operators, that violates barionic and leptonic numbers [27]. 

2 331 models with three identical fermion multiplets are usually call one family models. 
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2.1 Matching Conditions 

The general expression for the electromagnetic charge operator will be a linear combination of diagonal generators 
for the gauge group 331: 

Q = T 3 + Y = T 3 + 4=Ws + X (2.6) 
v3 

with Tj the Gell-Mann matrices normalized as Tr(TiTj) — The X operator for the abelian group U(l)x is 

proportional to the identity matrix 3x3. The hypercharge will be given by 

Y = -^=bT s + X. (2.7) 
V3 

6 is a known parameter that determines the class of 331 models to be considered [30 . 

The renormalization group analysis compares the couplings for different gauge groups at given energy scales, 
and models with symmetry breakings need relations for couplings at different energy regions which are called the 
matching conditions; they are extracted from the way in which the unbroken group is embedded into a bigger 
broken group. Also, in order to have all couplings in the same ground, all the generators should be normalized in 
the same way, and well normalized couplings are those which will converge in an unification point. 

Calling Y the well-normalized hypercharge operator, it will be proportional to the original one, 

Y = oY, (2.8) 

where a is the normalizing parameter so that the convergence of the running coupling constants at certain scale 
Mu is guaranteed. In the same way the operator X has a well-normalized X, i.e. X = cX, and the normalizing 
parameter for it, is given by Eq. (|2.7p requiring the same normalization for Y, T% which satisfies the following 
relation 

a 2 = \b 2 + c 2 . (2.9) 



Therefore, the parameter a should be such that 



a 2 > h 2 . (2.10) 




Then, the well-normalized hypercharge operator can be written as a function of the unknown parameter a, as 

(2.11) 

And from this equation, we obtain the following matching condition for the corresponding couplings |33j : 

a^ = ^l+(a 2 -^ax- (2-12) 

where ay, ax , a 3 L are related with U(\) Y , U{l) x and SU(3)l, respectively. 
The following relations must also be satisfied 

2 , 2 46 2 . 

a Y — a ay, a x = (a ^~' aXl 

u s = a 3C , a 2 L = a 3L - (2-13) 

where ax, ay and aih are related with U (l) x , U(l)y and SU(2)l, respectively. The third relation corresponds 
to the strong interaction where a s is associated with the Standard Model and a 3 c is related with the color part in 
the 331 model. Finally, the last relation corresponds to the embedding of SU(2)l into SU(3)l- 
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2.2 RGE analysis 

By replacing the relations described by Eqs. (|2.12[ [2~T3]) into Eq. (|2.1| . we can write the evolution for the RCC 
from the Z boson-pole Mz passing through a 331 symmetry breaking scale Mx, up to a certain Scale of unification 

M v , 



1 



4b 2 
3 



46 2 

oleu[M z ) 1 -—a 2 L{M z y 1 



by - ^fb 2L 



OL2l(Mz) 




(2.14) 

(2.15) 
(2.16) 



where the coefficients b s , b 2 L and by are related with SU(3)c, SU(2)l, U(l)y, respectively and they are calculated 
at energies in the range Mz < £t < Mx- The coefficients 63c, &3L, and bx, are related with SU(3)c,SU(S)l, 
and U(l)x respectively; they are calculated for energies in the range Mx < /i < Mjj- For our study, we need bi 
coefficients for energy scales below (and above) the symmetry breaking scale Mx, which will be given by SM (and 
331) degrees of freedom. Then for different models, we have different 6,'s for the intervals of energy scales which 
can change the running of the coupling constants. The input parameters from precision measurements are |31j 



Uem( m z) 



127.934 ±0.027, 



sin 2 6 W (M Z ) = 0.23113 ± 0.00015, 



a s (M z ) 

y 2L I 



a~l{M z ) 



0.1172 ± 0.0020, 
29.56938 ±0.00068. 



(2.17) 



The Mjj scale, where all the well-normalized couplings have the same value, can be calculated from (|2.15j) and 
(|2.16p as a function of the symmetry breaking scale Mx 



(M x 
= \~M~z 



exp < 2tt 



a s {Mz)- x - a 2L (Mz)- 1 



he - b. 



■3L 



(2.18) 



The hierarchy condition Mx < Mjj < Mpi an ck, must be satisfied. We shall however impose a stronger condition 
of Mjj < 10 17 GeV, in order to avoid gravitational effects. Hence, the hierarchy condition becomes 



M x < M v < 10 17 GeV 



(2.19) 



Such condition can establish an allowed range for the symmetry breaking scale Mx in order to obtain grand 
unification for a given normalizing parameter a. 

With a similar procedure, the expression for a 2 is found, and is given by 



45 2 
"3 



46^ 

T 



b — h 

^- + { aEMiMz)- 1 - ^ra^Mz)- 1 - Y 3 2L 

a s {M z )- 1 - a 2L {MzY 



ln 



( M x 
\M Z 



1 b,-b 



2L 



hi 



27T 63C - b?,L 

b 2 L 



( M x 
\Mz 



63c - b. 



■3L 



Oi 2 L {M z 



2tt 



(M x \ 




' 1 b s - b 2L 


\Mz) 


+ b 3L 


27T 6 3C - b 3L 



In 



( M x 
\M Z 



as(Mz)- 1 - a 2L (M z y 



b?,c — b 3 L 



> 



(2.20) 



In order to analyze the possibility of having unification at certain scale Mu with 331 as the only intermediate 
gauge group between the SM and Mjj scales, we could distinguish three scenarios of unification pattern 
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O The scenario with b 3 c b 3 u and (b 3 c — ^3l) 7^ (b s — b 2 u)- We shall call it the first unification pattern 
(1UP). In that case, we obtain an allowed region for the scale M x . It is carried out by combining the results 
of Section 12.21 The procedure to get the allowed interval for Mx is described in detail in Section 13.11 for the 
so called model A. 

© The scenario with 63c = b^u- We call it (2UP). In such a case Eq. (|2.18[) is not valid anymore, and we 
should go back to Eqs. (|2.15[[2~TB|) . For an arbitrary value of Mx the couplings a%c and «3L go parallel each 
other for energies larger than Mx- Therefore, the only possible way to still obtain unification is by setting 
&3c — ctsL at the scale Mx such that both couplings go together for scales above Mx- Unification with the 
third coupling could occur at any scale bigger than Mx- By equating Eqs (|2.15l I2.16| . we find a single value 
of Mx that makes the couplings a^c and a^u to converge. This convergence occurs at the scale 



M x = M z exp 



2ir [a 2L (M z )- 1 -a s (Af z )- 1 ] 



(b 2L - b s ) 



(2.21) 



It worths emphasizing that this scenario leads to a unique value of Mx and not to an allowed range. Finally, 
Eq. (|2.20p for a 2 must also be recalculated to find 



a 2 



F 2 (M x ) - H ln ( 



Mu \ 3 



„ s ,,,.-1 4fe 2 j b Y ~^-b 2L (2ir[a 2L (M z )- 1 -a s (M z )- 1 ] \ 
F l{ Mx) = a EM (M z ) i--a 2L (M z ) |_L _ lj 

F(M\ „ , M H b2L i ^[a 2L (M z )-^a s (M z )^] \ 

F 2 (M X ) = * 2L (M Z ) - — j — — j (2.22) 

© The 3UP with (fe 3 c — &3l) = (b s — b 2 i) 7^ 0; according to Eq. (|2. 18[) . the unification scale Mu becomes 
independent of Mx ■ 

The case (63c — &3l) = (b s — &2l) = 0, does not lead to unification as can be seen by trying to equate Eqs. 
(|2.15p and (|2.16[) . Since the first scenario is the most commom one, we shall only indicate when the other two 
scenarios appear. We will study non-SUSY and SUSY versions of the 331 models. In the case of SUSY models we 
shall consider two scenarios for the SUSY breaking pattern 

1. The SUSY Breaking Scenario at the Z— pole (ZSBS), in which the SUSY breaking scale is taken as Asusy ~ 
M z ■ Although this is not a very realistic scenario, numerical results do not change significantly with respect 
to the more realistic scenario with Asusy lying at some few TeV's. 

2. The SUSY Breaking Scenario at the M x scale (XSBS), with A SU sy ~ M x i.e. SUSY breaking at the 331 
breaking scale. 



3 331 Models 

Analogously to the SM, fermions will transform as singlets or in the fundamental representation of SU(3)l, and 
gauge fields in the adjoint representation. Assignment of U(\)x quantum numbers should be done ensuring a 
model free of anomalies. There are models with three different families necessary to cancel out the anomalies, and 
there are models with only one family and the other two are a copy of the first. We will take into account six 
different versions of the 331 model for the analysis of the unification scheme. 

The minimal spectrum necessary for symmetry breaking and generation of masses is given by [32) 

& ~ (1,3*, -1/3), 

02 ~ (1,3*, -1/3), (3.1) 
3 ~ (1,3*, 2/3). 

Where the quantum numbers are associated with SU (3) c , SU (3) L , and U (l) x respectively. The first multiplet 
acquires a vacuum expectation value (VEV) at Mx scale, breaking the symmetry as 331 — > 321; the other two will 
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be decomposed as singlets and doublets of SU(2)l, which will break the 321 symmetry. The spectrum transforms 
like 



, [fasti ~ (1,2*, -1/2), 

4,2 ~* U ~ (i.i.o)- ( } 

, [fosM ~ (1,2*, 1/2), 

* \4 ~ (1,1,1). (3 - 3) 

This is applied for models with 6 = 1/2, with tiny changes that do not affect our analysis 43, 44 . For models with 
6 = 3/2 this spectrum should be varied in order to get a phenomenological mass spectrum. For this case we have 
the following transformation for the triplets [3H 133 13S] : 

77~(1,3*,0), 

p~(l,3*,l), (3.4) 
X~(l,3*,-1). 

the 331 symmetry is broken when x acquires a VEV in the third component. The remaining scalars will transform 
in the 321 symmetry as 

^ hsM ~ (M*, -1/2), 
K ~ (1,1,1)- 
Jp SM ~ (1,2*, 1/2), 
P V++ ~ (M,2). (3 ' 6) 

The Model that we will call Model E needs also another scalar boson transforming under 6 representation of 
SU(3)l- Such an structure is required in order to give masses to the neutrinos [21)1 [341 135] . 

S ~ (1,6,0) (3.7) 

which has the following representation in the 321 symmetry: 

03SM ~ (1, 3, 1), 

cfosu ~ (1,2,1/2), (3.8) 
4>— ~ (1,1,-2). 

3.1 Model A 

The simplest anomaly- free structure for this gauge group is generated by taking b = 1/2; and the following spectrum 
is obtained [37] 

= (e-,i/ e ,JV?)£~(l,3*,-l/3), 

^2L = (^iV^ )!- (1,3*, -1/3), 

= (iV 4 °,E+, e +)l~(l,3*,2/3), 

Q L = («,<*,/>)£ ~ (3,3,0), (3.9) 



(3*, 1,-2/3), d£~(3*,l,l/3), D£~(3*,l,l/3). 



The other two families arc copies of the first one, and each family is free of anomalies. This particular 331 model 
could be embbeded in the E 6 gauge theory, but we shall assume that not necessarily such embedding occurs. 

From the general expression for coefficients bi Eq. (]2.2p . and using the quantum numbers assigned for each 
representation in this model, they can be expressed at energies below Mx in the following form 

by = y N * + 1 Nh + 1 £ y2(s) ' 

sing — s 

4 1 22 

b 2L = ^N g + -N H 



3 y 6 3 
4 

3 



b s = i;N g -ll. (3.10) 
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Figure 1: Allowed interval for Mx in model A. The solid line represents Mjj as a function of Mx- The horizontal 
short dashed line is the graphics for Mjj = 10 17 GeV. The long dash line represents the function Mjj = Mx- 
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a 2 
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MX 

Figure 2: Normalization parameter a 2 as function of the symmetry breaking scale Mx for model A. The intersection 
of this function with the constant function |6 2 gives a lower limit for Mx owing to the condition described by Eq. 



:oo 




100 100000. 1..10 8 1.-10 11 1.-10" 1.-10 17 

p(GeV) 



Figure 3: Running Coupling Constants for model A with Mx = 1-63 x 10 16 Gel^. The unification scale appears at 
M v ~ 10 17 GeV 
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where N g is the number of families, Njf the number of scalar doublets in SU(2)l, and the sum runs over scalar 
singlets. 

In the same way, the b^s for energies above the Mx are given by 



2 

:>> " ' 3' 
1 



b x = -N g + -, 



b 3L = 2Ng+--ll, 

2 

b 3C = 2JV fl -ll. (3.11) 
Therefore, from the spectrum showed in Eq. (|3.9p . we obtain 

(bY,b 2 L,b a ) = (^,-3,-7) (3.12) 



for energies below M x , and 



3 



'26 9 

(b x ,b 3L ,b 3C )=[-,--,-5) (3.13) 



for energies above Mx- 

First of all, by using the hierarchy condition Eq. (12.19j) . an allowed region for Mx can be found. In Fig. [T]we plot 
Mjj as a function of Mx, Eq. (|2.18|) . We also plot the line Mu = Mx, and the line Mu = 10 17 GeV, the intersection 
of these lines with the first plot, provides the allowed region given by the constraints Mx < My < 10 17 GeV 
according to Eq. (|2.19| . The intersection between the curve Mjj(Mx) from Eq. (|2.18p with Mu = 10 17 GeV gives 
the lower bound for Mx, and the intersection with Mjj = Mx gives the upper limit for Mx- Then from Fig. [1] 
and the hierarchy condition, it is obtained that 

1.63 x 10 16 GeU < M x < 2.05 x 10 16 GeU. (3.14) 

In Figure [2] we plot a 2 as a function of Mx from Eq. (|2.20|) . For this model 6=1/2 and from Eq. (12. 10[) the 
restriction a 2 > 1/3 is obtained. The horizontal line a 2 = 1/3 is also showed and the allowed values are above of 
this line. In this case the constraint for Mx is given by Mx < 2 x 10 14 GeV which is weaker than the previous 
restriction. 

From Eq. (|2.18p and the allowed interval for Mx of Eq. (|3.14p we get the region permitted for Mu i.e. 2 x 10 16 
GeV < Mu < 10 17 GeV. 

Now, for M x = 2.05 x 10 16 GeV (the highest allowed value) we obtain from Eq. (f2~T8|) that Mu « M x which 
is not a natural hierarchy. Instead, for Mx — 1-63 x 10 16 GeV (the lowest allowed value) the unification occurs 
at a more natural scale of 10 17 GeV. Then, we shall take the latter to plot the evolution of the running coupling 
constants. 

The evolution of the RCC described by Eqs. (|2.1[) . requires the results for the matching conditions Eqs. (|2.12[ 
I2.13p . Moreover, the following input parameters are necessary: ® The bt coefficients Eqs. (|3.12[ I3.13[) . ® The 
value of the 331 breaking scale Mx, and © the values of the couplings at Z-pole Eqs. (|2.17p . Assuming the lowest 
allowed value for the 331 breaking scale (Mx = 1-63 x 10 16 GeV) we plot the evolution of the running coupling 
constants in Fig. [3] The unification scale appears to be Mu = 10 17 GeV, in agreement with our previous results. 



3.1.1 SUSY Extension 

One possible SUSY extension in which the superpartners cancel the anomalies of their partners, arises by introducing 
new scalars which transform in the conjugate representations of the original ones. Under these considerations, and 
using Eqs. (|2.5p we obtain the supersymmetric bi coefficients 

(by,b 2L ,b s ) = (14,2,-3), (3.15) 
(bx,b 3L ,b 3C ) = (16,3,0). (3.16) 

We take first the ZSBS in which SUSY breaking is assumed to occur at EW scale d. Since in that case SUSY is 
not broken by going from EW to the unification scale, we should use the supersymmetric bi coefficients, at scales 



3 Our results are also valid for SUSY breaking at scales of the order of 1 TeV. 
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below and above the 331 breaking scale. Taking it into account, and using the hierarchy condition (|2.19p . as it was 
explained in the previous section, we get an allowed range for Mj, of 

1.27 x 10 8 GeV < M x < 2.76 x 10 13 GeV. (3.17) 

The bound coming from the a parameter Eq. (|2.10[) . yields 6.43x 10 6 GeV < Mx, which gives no further restriction. 
This allowed region for Mx leads to an allowed interval for Mjj obtained from Eq. (|2.18[) . getting 

2.76 x 10 13 GeV < M v < 10 17 GeV 

On the other hand, we shall consider the XSBS with SUSY breaking at the Mx scale. In that case, we should 
use the non-SUSY 6j coefficients for scales below Mx, and the SUSY bi coefficients for energies above Mx- From 
this we obtain the following allowed interval for Mx 

1.76 x 10 14 GeV < M x < 2.05 x 10 16 GeV, (3.18) 

for which the restriction a 2 > 1/3 is also satisfied. This allowed region provides a unification scale between 
2.05 x 10 16 GeU and 10 17 GeU which correspond to the values calculated from Eq. (|2.18p for the lower and upper 
values for Mx, respectively. Therefore, this is a possible scenario of unification with 331 as the unique intermediate 
gauge group, for both schemes of SUSY breaking. 

It worths noting that SUSY scenarios leads to much lower breaking scales for the 331 model, especially for low 
energy SUSY breaking. 

In another possible supersymmetric extension of model A, the scalar fields are not introduced explicitly. Instead, 
they appear as superpartners of the leptons, since they have the same representations and quantum numbers [45) . 
The supersymmetric bi coefficients read 

(by,b 2L ,b s ) = (10,0,-3), (3.19) 
(b x ,b 3L ,b 3C ) = (12,0,0). (3.20) 

For both schemes of SUSY breaking, the model correspond to the 2UP explained in Sec. (|2 . 2[) . so that Mx acquires 
a single value. The 331 scale appears at Mx — 10 21 GeV. It discards this SUSY version of model A, to get UCC 
under our scheme. 



3.2 Model B 

Again the parameter b takes the value 1/2, and in this case the 331 gauge theory is a family-symmetric model. For 
each family the spectrum is |38j : 

1>i L = {e-,v e ,Ei)Z~ (1,3*, -2/3), 

V^ 2 l = (Al°,£+X)I~(l,3*,l/3), 

V> 3L = {E-,N°,E-) T L ~ (1,3*, 2/3), 

e+ ~ (1,1,1), £+-(1,1,1), E+ L ~ (1,1,1), 

Ql = (u,d,I7)£~(3,3,l/3), (3.21) 

«i ~ (3*, 1,-2/3), d£~(3*,l,l/3), (3*, 1,-2/3). 

In this case, the 331 gauge group could be embbeded into the SU(6) U(l) x gauge theory, but we suppose it is 
not mandatory. This model has the same spectrum of model A at low energies. Therefore, it has the same values 
for the bi coefficients at low energies given by Eq. (|3.12l) 

(b Y ,b 2L ,b s )= (^,-3,-7). (3.22) 

For energies larger than M x , we have that and 63c are the same as in model A, since their values are ruled by 
the number of triplets, which coincide for both models. For bx, there is a difference coming from the additional 
singlets in the lepton sector. In this case, we obtain 

20 Ar 2 62 
b x = — N„ + - = — . 
3 9 3 3 



4 When a group of grand unification is present, protection from proton decay leads to a stronger bound for Mx by demanding 
Mxj > 2 X 10 16 GeV, in that case we find 1.27 X 10 8 GeV< M x < 1.4 x 10 9 GeV. 
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Hence, 

(bx,h L ,b 3C )= (f (3.23) 
Following the same procedure as for the Model A, we find the allowed range for Mx 

1.63 x 10 16 GeV < M x < 2.05 x 10 16 GeV. (3.24) 
In this interval, a 2 > 1/3 yields M x > 3.42 x 10 15 GeV, and M v lies in the interval [2.05 x 10 16 , 10 17 ]GeV. 

3.2.1 SUSY Extension 

As a first attempt, we take the model in Ref. [39], where scalars are considered as sleptons and two new superfields 
are introduced in order to give mass to the supersymmetric particles and cancel out the anomalies. They transform 
as 

</>4L ~ (1,3, 1/3), 4i ~ (1,3*, -1/3). (3.25) 
The spectrum is the same at low energies as SUSY model A, Eq. (|3.19|) . The supersymmetric b. L coefficients read 

(&r,&at,&a) = (10,0,-3), (3.26a) 
{bx,h L ,hc) = (24,3,0). (3.26b) 

For the ZSBS (Asusy ~ Mz), we are in the 3UP explained in Sec. (12. 2p so that Mu acquires a single value 

Mjj = 1.24 x l0 21 GeV. (3.27) 

which discards this SUSY version of model B with Asusy ~ Mz, to achieve unification with our breaking scheme. 
In the XSBS ( Asusy ~ Mx), the allowed region for the model is 

1.76 x lO u GeV < M x < 2.05 x 10 16 GeV. (3.28) 

the bound from a 2 is Mx > 1.44x 10 11 GeV giving no further restriction. Mu belongs to the interval [2.05 x 10 16 , 10 17 ] GeV. 
Hence, this SUSY version of model B with Asusy ~ Mx permits UCC under our assumptions. 

3.3 Model C 

In this model, the parameter b is also equal to 1/2, but the cancellation of anomalies is obtained with a number of 
families multiple of three [3D] . Its fermionic spectrum is 

il>l = {v a ,or,E-) T L ~ (1,3,-2/3), 

<4 ~ (1,1,1), (i.i.i), 

Qi = [d\u l ,U l ) T L ~ (3, 3*, 1/3), 

< c ~ (3*, 1,-2/3), (3*, 1,1/3), C/| c ~ (3*, 1,-2/3), 

Qi = {u 3 ,d 3 ,D)l~ (3,3,0), (3.29) 



l 3L 



(3*, 1,-2/3), ^-(3*, 1,1/3), D£~(3*,l,/3). 



where a — 1, 2, 3 label the three lepton families; while i = 1, 2 refer to only two of the quark families because the 
third one transforms in a different way. The spectrum is the same of the Model A at low energies, so that the bi 
coefficients are given by Eq. (I3.12p . For energies above Mx there are new contributions and the bi's yield 

22 



(b Y ,b 2L ,b s ) = ^-g-,-3,-7 

/42 13 \ 

(bx,b 3L ,b 3C ) = (y'-y'- 5 )- ( 3 - 3 °) 

Following the same procedure explained before, we find that the allowed region for Mx is 

2.05 x 10 16 GeU < M x < 3.15 x 10 16 GeU. (3.31) 

In this case a 2 > 1/3, leads to Mx < 9.57 x 10 18 GeV and provides an upper (and not a lower) limit because a 2 is 
a decreasing function of Mx- Of course, this constraint gives no further restrictions. Taking values for Mx from 
the interval in Eq. (pPTj) . the permitted interval for M v is [2.05 x 10 16 GeV, 10 17 GeV] . 
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3.3.1 SUSY Extension 

The extension is done adding new scalars transforming in the conjugate representation in order to cancel the 
anomalies. With this spectrum we obtain 

{bvMM = (14,2,-3) 
(bx,b 3L ,b 3C ) = (24,0,0). (3.32) 

At low energy scales, the coefficients coincide with the analogous in model A, Eq. (|3.15[) . For both schemes of 
SUSY breaking we lie in the 2UP getting a unique M x from Eq. (|2T2Tj) . For the ZSBS, M x becomes 

M x = 2.76 x 10 13 GeU. (3.33) 

If we fit Mu = 2.76 x 10 13 (=lx 10 17 ) GeV, we need a 2 = 3.26 (= 1.78) to obtain unification of the three couplings 
according to Eq. (j2~2"2")l . In general, for M v lying in the interval [2.76 x 10 13 GeV, 10 17 GeV] , we find that the 
a 2 parameter lies in the interval [1.78, 3.26]; and the condition a 2 > (4/3) b 2 is satisfied. Therefore, this model 
gives a consistent unification with 331 as the only intermediate gauge group. 

For the XSBS, M x = 2.05 x 10 16 GeV. From M v e [2.05 x 10 16 , 1 x 10 17 ] GeV, we get a 2 e [1.84, 1.97] and 
a 2 > (4/3) b 2 is fulfilled in all cases and the model is viable as well. 



3.4 Model D 

This model differs from model C just in its spectrum [411 142] . which is given by 







£~ (1,3,-1/3), 




<*£ 


~ (1,1,-1), 






Ql 


= {d\-u\D l ) 


£~(3,3*,0), 




< 


~ (3,1, +2/3), 


d£~ (3,1,-1/3), 


^ (3,1,-1/3), 


Ql 


= (d.3,u 3 ,U)l 


-(3,3,1/3), 




U C 3L 


~ (3,1,2/3), 


d c 3L ~ (3, 1,-1/3), 


VI ~ (3,1,2/3). 



(3.34) 

with ot = 1, 2, 3, running over the three families while i = 1, 2 runs over two families. The model differs from the 
others above in its spectrum at high energies. Therefore for the 331 gauge group the bi functions coincide with the 
values in the Eq. (I3.12| . 

(b Y ,b 2L ,b s )= ^,-3,-7j. (3.35) 
For energies larger than the symmetry breaking scale Mx we have 

(bx,bs L ,b 3 c)=(j-Y-^)- (3-36) 

These values give the following allowed region 

2.05 x 10 16 GeV < M x < 3.15 x W ie GeV, (3.37) 

and a 2 > 1/3 predicts Mx < 1-54 x 10 20 GeV, yielding no further restriction. The scale of grand unification lies in 
the interval 2.05 x 10 16 GeV< M v < 10 17 GeV. 



3.4.1 SUSY Extension 

Treating the model in the same way as model C, the bi coefficients at low energies are given in Eq. (|3.15[) , and for 
high energy scales they are given by 

(bx,hL,b 3C ) = (16,0,0). (3.38) 

Both schemes of SUSY breaking are in the 2UP. For the ZSBS wc find M x = 2.76 x 10 13 GeU. Taking M v € 
[2.76 x 10 13 , 1 x 10 17 ] GeV we obtain a 2 € [2.27, 3.26], and the basic restriction on a 2 is always accomplished. 
With XSBS, it is found that M x = 2.05 x 10 16 , and we get M v E [2 x 10 16 , 1 x 10 17 ] GeV for a 2 = [1.88, 1.97], 
respectively. 

Thus the SUSY version provides a possible scenario of unification for both schemes of SUSY breaking. 
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3.5 Model E 

In this model 6 = 3/2 and the cancellation of anomalies is obtained by the interplay of three families 26, 3~4l 135] . 
Its fermionic spectrum is 

r L = (e a ,v a ,e ca )T ~(1,3*,0), 

ql = {u\d\f) T L ~ (3,3,-1/3), 

qi = (d\u\ S f L ~(3,y,2/3), 

u c L a ~ (3,1,-2/3), -(3,1,1/3), si -(3, 1,-5/3), 

it ~ (3,1,4/3). (3.39) 

where a = 1,2,3 labels the families, and i = 2, 3 are related with two of them. 

Taking into account the whole possible scalar spectrum at low energy scales, including the scalars in the 6- 
dimensional representation, we can write the bi coefficients as 

b Y = ^N g + ±N H + ±3Y 2 (tmplet) + ± £ F 2 ( S ), 

singl — s 

4 11 22 
b 2L = -N g + -N H + -T R (tmplet)~ — , 



b s = ■zNg-U. (3.40) 



(bx,b 3L ,b 3C ) = (22,-^,-5). (3.42) 



where Y [triplet) and Tn(triplet) mean the quantum numbers of the SU(2)l scalar triplet in the 6-dimensional 
representation. Therefore the values of the bi coefficients read 

(b Y ,b 2L ,b s ) = (H.-H.-t), (3.41) 

17 

y 

With these bi coefficients we only find a lower limit from the hierarchy condition, 

6.87 x 10 13 Gey < M x (3.43) 

An upper limit is obtained from Eq. (|2.10p with 3 < a 2 , which is satisfied for Mx < 6.03 x l0 12 GeV, but both 
constraints are inconsistent and rule out the model under our assumptions. 

3.5.1 SUSY Extension 

The supersymmetric model is constructed adding scalar fields in the complex representation in order to avoid chiral 
anomalies from their superpartners. With this spectrum we obtain the following values for the bi coefficients 

(b Y ,b 2 L,b s ) = (28,6,-3), 
(bx,b 3L ,b 3C ) = (42,5,0). (3.44) 

We also consider an extended spectrum for the scalar sector. With these values we get only an upper limit for Mx 
in the ZSBS 

M x < 2.18 x lO s GeV. (3.45) 

But the restriction in the normalization parameter a in Eq. (|2.10p . requires Mx > 8.1 X 10 8 . This situation 
discards the model with only one intermediate breaking. 

For the XSBS, the hierarchy condition gives only a lower limit 6.87 x 10 13 GeV< Mx- But a 2 > 3 gives 
M x > 1.35 x 10 22 GeV, discarding the model. 
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3.6 Model F 

It is also possible to obtain cancellation of anomalies by a tiny variation of the fermionic spectrum from the model 
D of Pleitez-Tonasse |36j . This new spectrum will be given by 

r L = (v a ,e a -,E a +)f ~ (1,3,0), 

e a L + ~ (1,1,-1), (1,1,1). (3.46) 

where an exotic lepton has been introduced in the third entry of the triplet instead of the electron right-handed 
part, such that there are also leptonic singlets of SU(2>)l- 

The hi coefficients at low energy scales differ from those for the model E in the contribution due to the sextuplet. 
Therefore the values are 

(b Y ,b 2 L,b s ) = (26/3,-3,-7), 



{bx,b3L,b 3C ) = 




(3.47) 



With these coefficients, the allowed region for Mx is given by 

2.05 x 10 16 GeU < M x < 3.15 x 10 16 GeV. (3.48) 

The restriction from Eq. (|2.10p concerning the parameter a 2 , gives an upper bound of Mx < 4.1 x 10 17 GeV, i.e. 
no additional restriction. 



3.6.1 SUSY Extension 

The extension is done analogously to the model E, obtaining the following values for the bi coefficients 

(b Y ,b 2L ,b s ) = (22,2,-3), 
{bxMLMc) = (48,0,0). (3.49) 

The SUSY version is in the 2UP for both breaking schemes. In the ZSBS, Mx reads 

M x = 2.76 x 10 13 GeU. (3.50) 

In this case, the condition a 2 > 1/3 imposes a strong bound on the Mjj scale. From 

2.76 x 10 13 GeV < M v < 1.13 x 10 15 GeV 

we get a 2 = [1/3, 1.67]. Larger values of My leads to non allowed values of a 2 (less than 1/3). 

In the XSBS we have M x = 2.05 x 10 16 GeU. Using M v e [2 x 10 16 , 10 17 ] GeV we get a 2 £ [1.55, 1.82]. 
They satisfy the condition a 2 > 1/3. Hence, SUSY versions give a possible scenario of unification in our scheme. 



4 Some possible scenarios with groups of grand unification 

One of our main assumptions was that there is no necessarily a group of grand unification at the scale in which the 
coupling constants converge. In this way, the a 2 parameter becomes free, and indeed we could reverse the problem 
in a certain way, since we have in many cases an allowed region for a 2 and we could ask what groups of grand 
unification (if any) could lead to values of a 2 admitted by our scheme. Let us elaborate about this possibility 

After working some 331 models with one or three families, we can see that in some of them is possible to find 
a scale Mx and a normalization factor a 2 that gives unification of the coupling constants. In the absence of a 
group of grand unification (GUT) a 2 is not fixed by the group structure. Notwithstanding, UCC imposes some 
restrictions on a 2 , and by using these values we can look for simple groups containing a 331 group and fixing an 
a 2 belonging to the allowed interval mentioned above. Some good examples of GUT are E§, and SU (7). 

The group E 6 can be broken into SU (3) c <g> SU (3) L <g> SU (3) R or SU (6) L <g> SU (2). Models A and B shown 
in this paper could be embedded in this scheme. On the other hand, the 331 models of three families can be 
embedded into SU (7), thus, we shall study the latter scenario in more detail. 
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There are different versions of the SU(7) GUT that lead to 331 models according to the irreducible representa- 
tions (irrep) that cancel anomalies, and the definition of the electric charge or equivalently the linear combination 
of diagonal generators that define the hypercharge (i.e. the normalizing factor a 2 ). The combination of irre- 
ducible representations free of anomalies, permits to accomodate those 331 models with three families by the 
branching rules. The additional singlets or exotic fields acquire masses over the Mx scale, and do not affect the 
RGE. There are different combinations of anomaly free irreps. Taking into account the irreps , 5" 1 - ) [7], \t^?[21] 

and 5'™ 2 ^ 7 [35], where the subindex means the anomaly coefficient, the bracket coefficient means the dimension 
and the labels a, (3, 7 = 1, • • • ,7. Models can be classified if no irrep appears more than once i.e. a form like 
^(-i)q © ^7i) ® *(-2)q/37 i but in general the same irrep can be repeated such as 5 x ^(-i)q. © ^" 3 ^ © ^(2) 7 - ^ n 
the other hand, the branching rules according to (SU(3) C , SU(2)l) are given by 

W°$*^®$ a ^ = [(3*,l) + (l,2) + (l,l) + (l,l)] 7 e[(3*,l) + (3,2) 

+ (3, 1) + (3, 1) + (1, 1) + (1, 2) + (1, 2) + (1, 1)] 21 

© [(1, 1) + (3, 2) + (3, 1) + (3, 1) + (3*, 1) + (3*, 2) 

+ (3*,2) + (3*,l) + (l,l) + (l,l) + (l,2)] 35 (4.1) 

The electromagnetic charge of the particles can be chosen by defining the hypercharge which is a linear combination 
of the £7(1) factors of the SUJ7), and imposing the condition that the singlet of color have electromagnetic charges 
q = ±1/3, ±2/3, ±4/3, ±5/3(3 or that the electromagnetic charge of the leptonic singlets are ±1, 0. 

Assuming an unification scheme with the simple group SU(7) and passing through a 331 model with three 
families, by one step of symmetry breaking, we get the scheme 

SU{7) SU{3) C ® 517(3) L ® 17(1)* 

A ^ SU(3) C ® SU{2) L ®U{l) Y (4.2) 

the assignment of the electromagnetic charge is of the form 

Q = T 3L +aY a +[3Y b + 1 Y c = T 3L +aY 
a 2 = a 2 +/3 2 + 7 2 

where Y a , Y b , Y c have the same normalization as the T^l, Y generators, and they correspond to the U (I) 
abelian subgroups induced from different subalgebras of SU(7) D SU (n) L <g> SU (m) c ® U(l) a , where SU (n) L D 

SU(2)l <8 U(l) b and SU (m) c D SU(3) C ® U(\) c . If the fundamental representation is decomposed as in Eq. (|4.1|) . 
the most general assignment is 

Q = diag(q, q, q, b, a, a — 1, 1 — 3q — 2a — b). (4-3) 

And we can identify the charge of (3, 1) to be —1/3 and those of (1, 2) as (1, 0), then 

Q = diag(-l/3, -1/3, -1/3, b, 1, 0, -b) (4.4) 

where 6 = 0, ±1 if we have singlet leptons or b = ±1/3, ±2/3, ±4/3, ±5/3 if we have singlet quarks. 
From the following possible SU (7) maximal subalgebras 

Model I SU (4) c <g> SU (3) L <g> U (l) a 
Model II SU (3) c O SU (4) £ O U (l) a 
Model III SU(6) L ® U(l) a -> SU (3) c <E> SU (3) L <g> U (l) c ® U (l) a 

we can settle three different assignments for the hypercharge according to the scheme in Eqs. (|4.3l 14. 4|) . For 
example, Model I is known in the literature as the Pati Salam model. The generators are defined by 

Y a = ^diag (1/4, 1/4, 1/4, 1/4, -1/3, -1/3, -1/3) 
Y b = J\diag (0,0, 0,0, 1/2, 1/2,-1) 

Y c = \ \diag (1/3, 1/3,1/3,-1,0,0,0) (4.5) 



5 The values 4/3, 5/3 correspond to exotic quark charges that arise in some 331 models, like model E. 
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Table 1: Hypercharge definition using the embbeding of £7(1) into subalgebras of SU(7) for models I, II, III 

In this case Y b corresponds to the generator Tg of SU(3)l', and a linear combination of Y a , Y c corresponds to 
U(l)x of the 331 models. For model II the generators are defined by 



Y a = \j-diag (1/3, 1/3, 1/3, -1/4, -1/4, -1/4, -1/4) 



y° = 

yc = 



diag (0,0,0,0,1/2,1/2,-1) 
diag (0,0,0,-1,1/3,1/3,1/3) 



(4.6) 



And for model III the generators are defined by 



ya _ 

Y b = 

yc _ 



-dio fl (1/6, 1/6, 1/6, -1, 1/6, 1/6, 1/6) 

r 

12 
3"' 



diag (1,1, 1,0, -1,-1,-1) 
diag (0,0, 0,0, 1/2, 1/2,-1) 



(4.7) 



By choosing the values for the free parameter b = 0,±1 for leptons, and b = ±1/3, ±2/3, ±4/3, ±5/3 for quarks, 
the coefficients a, f3, 7 that determine the hypercharge can be obtained, and they are displayed in table [TJ On the 
other hand, table [2] shows the a 2 normalization factor of the hypercharge according to the b factor. 

The next step is to check whether the three family models consider here (models C, D, E, F) can be accomodated 
properly in a SU (7) GUT. It can be done by comparing the allowed region for a 2 in tabled with the values of a 2 
in table [2] The second column shows the allowed interval of a 2 obtained in previous sections in the absence of a 
grand unification group, the third column shows the allowed interval by taking into account additional restrictions 
coming from proton decay (Mu > 2 x 10 16 GeV), that arise when a GUT theory is introducecH. 

The restrictions displayed in table[3]for non SUSY versions of models C, D, F; show that they can be embedded 
in SU (7) since the bound a 2 > 1/3 is accomplished by all the values of a 2 in table [2] Model E (the only one with 
exotic charges studied here) has been ruled out from phenomenological grounds under our scheme. On the other 
hand, from tabled we see that the a 2 parameter is strongly restricted in the SUSY versions of these models, and 
not all of them can be accomodated in SU (7) according to tabled 
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b 







^ ~ 1.67 


±1 


¥ ~ 5.67 




^2.11 




4f — 3. 44 


-p- 


^ ~ 8.78 




^ ~ 12.78 



Table 2: Hypercharge definition using the embbeding of Z7(l) into subalgebras of SU(7) for models I, II, III 





Allowed interval for a 2 with no 
restriction from proton decay 


Allowed interval for a 2 with 
restriction from proton decay 


non SUSY C, D, F 


> 1/3 


> 1/3 


C (SUSY ZSBS) 


[1.78, 3.26] 


[2.07, 3.26] 


C (SUSY XSBS) 


[1.84, 1.97] 


[1.84, 1.97] 


D (SUSY ZSBS) 


[2.27, 3.26] 


[2.47, 3.26] 


D (SUSY XSBS) 


[1.88, 1.97] 


[1.88, 1.97] 


F (SUSY ZSBS) 


[1/3, 1.67] 


Excluded 


F (SUSY XSBS) 


[1.55, 1.82] 


[1.55, 1.82] 



Table 3: Allowed interval for the normalization a 2 for models with three families. In the second column, no 
restrictions from proton decay are taken into account, while in the third column they are. Models E (SUSY and 
non SUSY) were discarded from phenomenological grounds. 

5 Conclusions 

We have studied the possibility of having a 331 model as the only intermediate group between the electroweak 
scale and the scale Mjj in which the three coupling constants unify. We assume that there is no necessarily a 
simple gauge group at the M\j scale. From the analysis of the renormalization group equations (RGE) we examine 
different 331 models of one and three families as well as their supersymmetric extensions. 

Specifically, we are supposing that the three different couplings unify at certain scale Mjj, and a symmetry 
breaking to a gauge group SU {2>)l®U (1)x occurs at a lower scale Mx- Then, a second breaking occurs at the Mz 
scale to arrive at the SM gauge group. We are also assuming that all particles beyond the SM are getting masses 
of the order of Mx- Other analyses can be done by supposing that the particle contents are in different thresholds. 
For our analysis of the RGE, the following conditions are taking into account: O The hierarchy condition given by 
Eq. (f2T9f i.e. M x < M v < 1 x 10 17 GeV, should be satisfied. © The condition described by Eq. (|2~10| should be 

6 There is another correction coming from the new spectrum introduced, but we have assumed that the new heavy modes do not 
alter the RGE significantly. 





Allowed interval for Mx 


Model 


No SUSY 


NSBS 


SBS 


A 


[1.63 x 10 lfi , 2.05 x 10 16 ] 


[1.27 x 10 8 , 2.76 x 10 13 ] 


[1.76 x 10 14 , 2.05 x 10 16 ] 


B 


[1.63 x 10 16 , 2.05 x 10 lfi ] 


No allowed 


[1.76 x 10 14 , 2.05 x 10 16 ] 


C 


[2.05 x 10 16 , 3.15 x 10 16 ] 


= 2.76 x 10 13 


= 2.05 x 10 16 


D 


[2.05 x 10 16 , 3.15 x 10 16 ] 


= 2.76 x 10 13 


= 2.05 x 10 16 


E 


No allowed 


No allowed 


No allowed 


F 


[2.05 x 10 16 , 3.15 x 10 16 ] 


= 2.76 x 10 13 


= 2.05 x 10 16 



Table 4: Summary of the results obtained in the paper for the six different 331 models studied, as well as their 
corresponding SUSY extensions 
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fulfilled, where a is the normalization parameter of the hypercharge Y, that leads to the unification of the couplings 
at certain scale. © There is no neccesarily a grand unification gauge group at the Mjj scale. 

In the case of the supersymmetric extensions, the RGE analysis depends on the specific supersymmetric version 
but also on the supersymmetry breaking scale. In particular, we assume two possible SUSY breaking scenarios: 
SUSY breaking at electroweak scale and SUSY breaking at the scale of 331 breaking. Although SUSY breaking at 
electroweak scale is not a realistic scenario, numerical results do not change significantly with respect to the more 
realistic framework with SUSY breaking at some few TeV's. 

Based on the criteria explained above, we can either find an allowed interval for the 331 breaking scale Mx 
or rule out the model as a possible grand unified theory under the scheme described above. We summarize the 
results found in this paper in table [4] for the six different 331 models studied here as well as their SUSY extensions 
(some SUSY extensions that are ruled out have not been included). In many cases, we see that SUSY versions tend 
to give lower allowed values for the 331 breaking scales than their non SUSY counterparts, it could make SUSY 
extensions easier to test from the phenomenological point of view. It worths emphasizing that taking into account 
that the spectrum out of the SM lies at the Mx scale, when assuming that the SUSY breaking occurs at the Mx 
scale we are in a natural scenario for split Supersymmetry [46j-|47j. 

On the other hand, we see that models C and D predicts the same allowed intervals for Mx in SUSY and non 
SUSY versions. However, SUSY versions predict different values for the normalization parameter a 2 . 

In addition, by finding the allowed values for a 2 we can proceed to see what groups of grand unification could 
give a value of a 2 lying in the allowed range. Although the introduction of a group at Mjj scale introduces new 
singlets that could lead to tiny changes in the RGE analysis from Mx to Mjj, this procedure helps us to figure 
out what scenarios of grand unification (if any) are possible under our scheme. In particular, we find that some 
331 models with three families, can be properly embedded in a grand unification scenario with SU (7), especially 
in non susy frameworks. 

As a matter of perspectives, two loops analysis for RGE can be carried out [JS] m order to fit the Mjj scale better, 
especially taking into account the uncertainty in the starting coupling constants at electroweak scale. However, 
significant changes in the allowed regions obtained here are not expected for a well behaved perturbative regime. 
On the other hand, it is possible to unify the coupling constant associated to U (1) at string scale instead of the 
M v scale @9]-[50]. 

Finally, it worths saying that the non allowed models are not neccesarily ruled out. We only can say that they 
cannot produce UCC under the scheme in which 331 is the only gauge symmetry between the Mjj and SM scales. 
For instance, it could be possible that they achieve unification by either introducing new physics, or requiring extra 
breaking steps from the Mjj scale to SM scale. By introducing new physics, it is possible to get a lower Mx scale 
of the order of TeV accesible to the LHC. 
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